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We apply a hybrid approach which combines loop and Fock quantizations to fully 
quantize the linearly polarized Gowdy T 3 model in the presence of a massless scalar 
field with the same symmetries as the metric. Like in the absence of matter content, 
the application of loop techniques leads to a quantum resolution of the classical 
cosmological singularity. Most importantly, thanks to the inclusion of matter, the 
Q \ homogeneous sector of the model contains flat Friedmann- Robertson- Walker (FRW) 

solutions, which are not allowed in vacuo. Therefore, this model provides a simple 
setting to study at the quantum level interesting physical phenomena such as the 
effect of the anisotropies and inhomogeneities on flat FRW cosmologies. 



o 

(N 
O 



o 
cr. 



> 
(N 

m 

(N 

o 



X 



PACS numbers: 4.60.Pp, 04.60.Kz, 98.80.Qc 



* mercedes@aei.mpg.de 

t daniel. mart in@iem. cfmac .csic.es 

* mena@iem.cfmac.csic.es 



2 



I. INTRODUCTION 

Loop quantum cosmology (LQC) [H-|H is a quantization of cosmological models inspired 
in loop quantum gravity ideas and methods p, 0], in which the geometry has a discrete 
quantum nature. The first model successfully quantized to completion in LQC was the flat 
FRW model minimally coupled to a massless scalar field, whose dynamical analysis shows 
that a quantum bounce replaces the initial singularity 0] • The resolution of the cosmological 
singularity is a robust property of the theory [gl, [9| , owing to the polymeric representation 
adopted for the geometry, and it is also achieved in the rest of homogeneous models quantized 
so far in LQC (see for instance 1Q-18[ and references therein). 



In order to allow for the presence of inhomogeneities within the framework of LQC, re- 
cently a hybrid approach to the quantization has been dev elop ed in the example of the 



simpler case: the Gowdy T 3 model with linear polarization jl9l-|22|. This is a midisuper- 



space with three-torus spatial topology that contains inhomogeneities varying in a single 



direction 23 



The introduced hybrid approach combines the techniques of LQC with those of the Fock 
quantization for reduced models in which only global constraints remain to be imposed at the 
quantum level. The phase space is split in homogeneous and inhomogeneous sectors. The 
former is described by the degrees of freedom that parameterize the subset of homogeneous 
solutions, and the second one is formed by the rest of degrees of freedom. In the quantum 
theory, the inhomogeneous sector is represented a la Fock, in order to deal with the field 
complexity, while the homogeneous sector is represented following LQC, with the aim at 
obtaining a quantum model with no analog of the classical cosmological singularity. The 
approach assumes a hierarchy of quantum phenomena, so that the most relevant effects 
of the loop quantum geometry are those that affect the homogeneous degrees of freedom. 
In the case of the quantized Gowdy model, the homogeneous sector coincides with the 



phase space of the Bianchi I model, which has been extensively studied in LQC [14H16I. |24 
Concerning the inhomogeneous sector, the requirement that the conventional description 
for the inhomogeneities should be recovered when the quantum geometry effects of the 
homogeneous sector are negligible and that this description respect unitarity selects the Fock 



quantization of Refs. |25H27j without ambiguity. In fact, with the commented requirement, 



it has been shown that this is the unique satisfactory Fock quantization that the totally 



deparameterized Gowdy T 3 model admits [28J, |29 

Our aim is to further analyze inhomogeneous cosmologies in LQC by means of this hybrid 
quantization, now allowing for the presence of matter. In order to do this, we will include 
in the Gowdy T 3 model a minimally coupled massless scalar field with the same symmetries 
of the geometry. Choosing suitable field parameterizations for the inhomogeneities of both 
gravitational waves and matter, the corresponding field contributions appear in the con- 
straints in the same way [30|, 3l||. As a consequence, the uniqueness results of Refs. 28, 29 



for the Fock quantization apply to the nonvacuum case as well, and hence we have at our 
disposal a preferred Fock description also for the inhomogeneities of the matter field. 

The interest of this work lies not only in the fact that it provides a complete quantization 
of a cosmological model with an inhomogeneous matter field in the framework of LQC, but 
also in that it means a further step towards the quantum analysis of physical inhomogeneities 
in cosmology, in the sense that these inhomogeneities propagate on a geometry not very 
different from that of our universe. Indeed, thanks to the inclusion of matter, now the 
homogeneous sector of the model (nonvacuum Bianchi I) admits flat FRW cosmologies as 
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a subset of solutions, namely the isotropic ones, and it is widely known that the observed 
universe can be approximated at large scales by a spacetime of this type. Therefore, it is 
natural to compare the dynamics of our inhomogeneous model with that of the flat isotropic 
model, and analyze the quantum effects that anisotropies and inhomogeneities produce over 
a hypothetical FRW-like background. In particular, this setting would allow us to investigate 
questions like the robustness of the quantum bounce scenario of LQC when inhomogeneities 
are included, or modifications to the evolution of the matter inhomogeneities when quantum 
geometry effects are taken into account. 

Let us mention that, owing to the isometries of the Gowdy model, this family of space- 
times presents a particular subset of solutions with local rotational symmetry (LRS), in 
which the two scale factors of the directions of homogeneity coincide. Therefore, in order to 
simplify the analysis, it is convenient to focus on this kind of solutions, which we will call 
LRS-Gowdy cosmologies in the following. We will carry out this LRS reduction at the quan- 
tum level using an adaptation of the (so-called) projection procedure introduced in Ref. [16 
to pass from the loop quantized Bianchi I model to the loop quantized FRW model. 

This work is intended as a first contribution to the analysis of the Gowdy system with 
matter. Specifically, here we quantize to completion the model, putting special attention to 
the new features that the consideration of the matter field introduces in comparison with 
the vacuum case. We also present a general discussion of the lines of attack that can be 
pursued to extract the physics from our quantum model. We leave for a future work a more 
rigorous and deeper study of this underlying physics and its consequences. The structure of 
this paper is as follows. The classical model is described in Sec. [Til In Sec. IHH we carry 
out its quantization, promoting the constraints to operators and characterizing the physical 
Hilbert space. We also show there how to reduce the quantum model to the corresponding 
LRS-Gowdy counterpart. Finally, in the concluding Sec. HVJ we point out the possibilities 
that this model provides to analyze quantum phenomena in cosmology and to reach physical 
predictions. 



II. CLASSICAL MODEL 

The linearly polarized Gowdy T 3 cosmologies are globally hyperbolic spacetimes with 
three-torus spatial topology and two axial and hypersurface orthogonal Killing vector fields 



23J. We provide this system with a minimally coupled massless scalar field, $, with the same 
symmetries. We use global coordinates {t,9,a, 8}, where 9, a, 5 G S 1 , such that the Killing 
fields are d a ,ds. Then, all the fields (metric and matter) only depend on the coordinates t 



and 9. We reduce the system by performing a partial gauge fixing, as in Refs. [20|, |21|, [26 
As a result, the gravitational sector of the phase space turns out to be described by two 
pairs of canonically conjugate point-particle variables (they do not depend on 9) and by one 
field, together with its canonical momentum. We expand the fieldlike variables in Fourier 
series in the coordinate 9 and split the phase space into two sectors: one formed by all 
the homogeneous degrees of freedom (the two point-particle gravitational variables and the 
zero modes of both matter and gravitational fields, together with their momenta) and the 
other formed by the nonzero modes of the two fields of the system and of their conjugate 
momenta. We call them homogeneous and inhomogeneous sectors, respectively. 

In the totally deparameterized model, there is a particular field parameterization of the 
metric in which the gravitational wave content is described by a field which behaves exa ctly 



as the matter field $, namely, as a massless scalar field propagating in 2+1 gravity 32 . 
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Nonetheless, this description does not admit any Fock quantization with unitary dynamics 



29l . l33l . l34j . For that, it is necessary to apply a time dependent canonical transformation 



on both fields [26|, |29j. The resulting gravitational and matter fields, which we call £ and ip 
respectively, follow the equation of motion of a free scalar field with a time dependent mass in 
a static spacetime of 1+1 dimensions. Consistent with our restriction to the inhomogeneous 
sector, we consider both fields already devoid of zero modes. Now, introducing for these 
fields creation and annihilation-like variables defined like one would naturally do in the case 
of free massless scalar fields, one reaches a Fock quantization whose evolution is indeed 



unitary 25|, |26| and such that the vacuum is invariant under S 1 translations, which is the 
gauge group of the reduced system. Moreover, it has been shown that these two natural 
properties of unitary dynamics and vacuum invariance in fact pick up this Fock quantization 



as the unique acceptable one, up to unitary equivalence |28|, |29[. So, taking into account this 
result for the totally deparameterized model, we adopt the suitable field parameterization of 
Refs. M,M both for the gravitational and matter inhomogeneities of our current model, 



and describe them in terms of the creation and annihilation-like variables mentioned above, 
in order to eventually carry out the corresponding Fock quantization of the inhomogeneous 
sector. We will call these variables (flm'*, al^), with m G Z — {0} and a = £, ip. 

On the other hand, the homogeneous sector describes Bianchi I cosmologies with spatial 
three-torus topology and with a minimally coupled homogeneous massless scalar field, given 
by the zero mode of $. From now on, we call <fi this homogeneous matter field and 
its momentum. Since the homogeneous sector is to be quantized using LQC methods, 
we describe the gravitational variables of this sector in the Ashtekar-Barbero formalism. 
Using a diagonal gauge, the nontrivial components of the densitized triad are Pj/An 2 , with 



j = 6, a, 5, whereas those of the su{2) connection are Cj/2n (see e.g. [15]). These variables 
satisfy {ci,Pj} = SnG'jdij, where 7 is the Immirzi parameter and G is the Newton constant 
(throughout the text, we set the speed of light equal to the unity). 

Two global constraints still remain on this reduced system: the spatial average of the 
densitized Hamiltonian constraint, C, and the generator of S 1 translations, Cg. On the 
one hand, C can be split into two terms, C = Chom + C imi , the former involving only the 
homogeneous sector. Then, C in h couples the homogeneous gravitational sector with both 
the gravitational and matter inhomogeneous sectors through two identical terms, one per 
field, Cinh = Cf nh + C? nh , where Cf„ h denotes the corresponding coupling term for the vacuum 
case [2C|. The homogeneous term Chom is the densitized Hamiltonian constraint of the 
Bianchi I model minimally coupled to a homogeneous massless scalar field [la, E^] • Thanks 
to the presence of matter, the classical Bianchi I model admits solutions of the FRW type. On 
the other hand, as it happens to be the case with the inhomogeneous term C imi , Cg is the sum 
of two identical contributions, Cg = Cq+Cq, where Cq denotes the analog constraint in vacuo 
[20| . We see that, with our choice of variables, matter and gravitational inhomogeneities 
contribute in the same way to the constraints, and then it is straightforward to promote Cg 
and C to operators following the hybrid quantization developed for the vacuum case. 



III. QUANTUM MODEL 
A. Kinematics and constraint operators 



The quantization of the system starts with the introduction of a kinematical Hilbert 
space where the basic variables are represented as operators and where the constraints are 
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imposed quantum mechanically. For this kinematical Hilbert space, "Hkm; a natural selection 
is the tensor product of the kinematical Hilbert space of the gravitational sector "H|^ v times 
the kinematical Hilbert space of the matter sector "H™^ tt . Both of these spaces are in turn 
the tensor product of two spaces corresponding to homogeneous and inhomogeneous sectors, 
respectively. Physically, the nontriviality of the system comes from the couplings introduced 
at the moment of imposing the quantum constraints. 



For the gravitational sector, then, we carry out the hybrid quantization of Refs. [21], |22 
namely WE ™ is the tensor product of the kinematical Hilbert space of the Bianchi I model 



in LQC [16|, |22|, %^ n , times the standard Fock space for the inhomogeneities, J 7 ^, defined 
in terms of the annihilation and creation variables previously described for the field £. The 
homogeneous matter sector, on the other hand, is formed by the zero modes of the massless 
scalar field and its momentum, determined by <ft and P^. In analogy with the nonvacuum 
cases analyzed in homogeneous LQC (in particular the Bianchi I model minimally coupled 
to a massless scalar [lf|), we take the standard representation for these variables, choosing 
L 2 (IR, dtp) as the Hilbert space. Finally, since matter and gravitational inhomogeneities have 
identical behavior, the kinematical Hilbert space accounting for the matter inhomogeneities, 
J rip , is totally analogous to J-^. Summarizing, 

H kin = H^ n ®L 2 (M,#)®J^®J^. (3.1) 

For the inhomogeneous sector, the chosen representation is obtained by promoting the 
classical variables a"m* and dm\ with m G Z — {0} and a G {£, y?}, to creation and annihi- 
lation operators, a™ and dm , respectively. With them, it is straightforward to construct 
the quantum counterpart of the constraint Cg, for which we choose normal ordering. The 



result is 20, 21 



Co = £ mXi + £ mil, * m = o£)to&) - oWtfiW . (3.2) 



m=l m=l 



The same strategy is adopted when representing the inhomogeneous contributions to the 
coupling terms C? h . It turns out that the inhomogeneities of the field a (£ or tp) appear 
in Cj" h only via two different quadratic combinations, Hq and H^ t , whose normal ordered 



quantum counterparts are [20j, [21 



oo oo 1 

HS = miV m , = ~ (K + a&^SI + . (3-3) 

m=l m=l 

with = Om'^m + a_l^-m' The above operators X m , i^Q , and if^ t act nontrivially 
on T a and have as a common dense domain the space of n-particle states. We call the 
number of particles of the field a in the mode m. 

On the other hand, for the homogeneous sector, the basic matter variables are represented 
by the operators 0, which acts by multiplication, and = —iHd^, while for the gravitational 
part we adopt the operator representation discussed in detail in Ref. (HJ (see also Ref. (22 |). 
adhering to the improved dynamics scheme put forward by Ashtekar and Wilson-Ewing [16( 
(and which was called "case B" in Ref. (ll[). Let us briefly review this quantization scheme. 
First we recall that, on Ti^, the operators pi (i = 9, a, 5), which represent the nontrivial 
coefficients of the densitized triad of the Bianchi I model, have a discrete spectrum equal 
to the real line. The corresponding eigenstates, \pg,p a ,p$), form an orthonormal basis (in 
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the discrete norm) of H^L- Owing to this discreteness, there is no well-defined operator 
representing the connection, but rather its holonomies. The representation of the matrix 
elements of these holonomies incorporates the so-called improved dynamics prescription, 
which states that there exists a dynamical (state dependent) minimum length /2j for the 
straight edges in the zth-direction along which the holonomies are computed. We use the 
specific improved dynamics prescription put forward in Ref. [3]. Then, the elementary 
operators which represent the matrix elements of the holonomies, called Jx^, produce all 



a constant shift in the physical Bianchi I volume [16|, [21[. The resulting action of on 
the states \pe,Pa,Ps) is quite involved. In order to simplify the analysis, it is convenient to 
relabel the basis states in the form \v, A CT , \$), where v is an affine parameter proportional to 
the volume of the compact spatial section, such that any of the operators J\±p, i (i = 0, a, 5) 
causes a unit (positive or negative) shift on it. The parameters A, are all equally defined 
in terms of the corresponding parameters pi, and verify that v = 2\o\ a \$ (see the explicit 
definitions in Ref. |16||). 

Employing the basic homogeneous gravitational operators pi and Jx±p, i , we can complete 
the construction of the constraint operator C = Chom + Cinh exactly in the same way as in 
the vacuum case j2l[. This densitized Hamiltonian constraint operator is formed by 
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JJ mt ' 11 mt / ) 



(3.4) 



(3.5) 



with i,j G {8, 5, a}. Here, [1 / 1 jo^ | 1 / 4 ] is a regularized triad operator which has a diagonal 
action on the considered basis of states. On the other hand, the operator 0j is the quantum 
counterpart of the classical quantity CiPi and its action on the basis states is highly nontrivial. 
In particular, Oj and Qj do not commute for i ^ j. We will not give here the explicit action 
of these operators on our basis states (which can be found in Ref. (ll[). Instead, in the 
following section, we will write down explicitly the general equation that must be satisfied 
by the solutions of the quantum densitized Hamiltonian constraint. 

The above constraint operator leaves invariant certain subspaces of "Hkm, which provide 
superselection sectors 0, 2l[. When symmetrizing C, we have chosen a specific factor or- 



dering which leads to superselection sectors which are particularly simple and with most 
convenient properties. More precisely, instead of considering TL^, we can restrict the ho- 
mogenous gravitational sector to be the completion with respect to the discrete norm of 
the space spanned by the states \v, A CT , \$) such that v, \ a , and \$ belong to an octant, for 
instance v , X a , A# > (case on which we will focus our attention from now on), and with v 
belonging then to any semilattice C e of step four included in M + : 



£ e = {e + 4fc; keN}. 



(3.6) 



In this expression, e is any number in the interval (0,4], and provides the minimum value 
that v takes. In addition, given e, the labels A a (a = a or 5) are restricted to sectors of the 
form A a = X*uj t , where the A*'s are any two fixed positive numbers and ui e runs over the 
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following numerable and dense subset of 1R + : 

/ „ + 9™ , \ Pfc 1 

(3.7) 




Here rrik,nk,Pk G N, and zfZ when e > 2, while z = otherwise 21 . 

Once we restrict the study to any of the above superselection sectors, the null eigenspace 
of the homogeneous densitized triad operator (which is a proper subspace of "Hj^J ceases to 
be included in our theory. As a consequence, there is no analog of the classical cosmological 
singularity in the quantum model anymore. In this sense, it is ensured that the singularity 
is resolved, already at the kinematical level. 



B. Physical Hilbert space 

Once we have constructed the constraint operators, we can proceed to determine the phys- 
ical states, which must be annihilated by these constraints. Notice that the two constraint 
operators commute and can hence be imposed consistently. 

Let us consider first, e.g., the S 1 symmetry generated by Ce, which amounts to the 
following condition 



m=l 



m 



(At + X£) = 0, X «=n a m -n a _ mJ a = t,<p. (3.8) 



The states that satisfy this condition form a proper subspace of ® J 7 ^, which we call J-" p . 

The Hamiltonian constraint operator imposes a more complicated condition, mainly be- 
cause of the nontrivial actions of both ©j on the homogeneous gravitational sector and H^ t 
on the inhomoereneous sector [HQ- For our purposes here, it suffices to make explicit the 



action of the Hamiltonian constraint operator on just the homogeneous sector. With this 
aim, it proves convenient to introduce an alternate labeling of the basis states of "H^. The 
new labeling is given by \v, A, T), where A = In^X^Xs) and T = ln(Xs/X a ). Next, we expand 
a general state in this basis: 

|*)= \S>(v,A,r))®\v,A,r). (3.9) 

v,A,f 

Here, v, A, and T take values in the corresponding superselection sectors. Let us clarify that 
the kets \^(v, A, T)) are actually not wave function coefficients, but rather states inasmuch 
as we have not expanded | in a basis of the whole kinematical Hilbert space, but only of the 
homogeneous gravitational sector. On the other hand, based on our experience with other 
similar cosmological models, the states \^f(v,A, Y)) for the solutions of the Hamiltonian 
constraint are not expected to be normalizable in L 2 (R, dip) <S> J 7 ^® J 7 ^, but rather to belong 
to a larger space from which one should construct the physical Hilbert space of the theory. 
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Acting with C on and projecting over (v, A, T|, we obtain: 



|*(«,A,T)>+ E E + A, T)> 

«e{o,4} se{+,-} 

4/3e 2A 6 2 ( V )[^ t + ^n t ] E E fo2 ^ + 

KG{0,4} sG{+,-} 



+ 



/3e 



2A 



'Hl + HS] |*(t;,A,T)>=0. 



(3.10) 



Here, (3 = [Gh/(lfa 2 j?AW/ 3 , with A denoting the minimum nonzero eigenvalue allowed 
for the area in LQG [ly, |2l| , and we have defined 



b(v)= 

<0) = 



v- 1 



3 17Tk/4 



s2 + d | a/w fs/c + f | s 1 + sgn ( s 



2 + - 
2. 



+ V 



(3.11) 
(3.12) 



On the other hand, the objects \^ S K (sk + v, A, T)), are linear combinations of six contribu- 
tions in the form 

|**(s« + v,A,T)>= E (\*(sK + v,A + w v (s2),T + rw v (s2))) 

re{l,-l} 

+ |\P (sk + v, A + w v (sk), T + r w v (sk) — 2r w v (s2))) 

+ I^JJ (sk + d, A + w v (sk) — w v (s2), T + r w v (sk) — r u;„(s2)))J , (3.13) 

where w v (sn) = ln(sn — ln(t>). The two last lines of Eq. ( 13 . 1 j) correspond to the action 
produced by C in h, where we have introduced the notation 

\W%(sk + v,A,T)) = E (l* (sk + v,A + w v {sk),T + tw v (sk))) 

re{i-i} 

+ |* (sk + v,A + w v (sk),T + rw v (sn) -2rw v (s2)))\ (3.14) 



Condition (I3.10p . coming from the constraint, is a difference equation in the variable v and 
can be seen as an evolution equation in this variable. In the vacuum model, it has been 
proven that, formally, a(n infinite but countable) set of initial data on the section given by 
the minimum value of v, Vmip ~ e £ (0, 4], completely determines a solution of the densitized 
Hamiltonian constraint [2l[ 22| • Since the Hamiltonian constraint of our model and the one 
in vacuo have identical structure, the above result applies also to our case. This property 
allows us to identify the physical Hilbert space of the system, that we call "H p h y s, as the 
Hilbert space of these initial data. 

The resulting physical Hilbert space, taking into account condition (13. 8p as well, is given 
by "Hphys = %ph ys ® L 2 {R, d4>) ® J- v , where "Hph ys is the physical Hilbert space of the Bianchi 
I model determined in Ref . [22j . As discussed in that reference, the inner product that pro- 
vides this Hilbert space structure on the space of initial data is obtained by the requirement 
that the complex conjugation relations between a complete set of classical observables turn 
into adjoint relations between the corresponding operators. 
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C. Projection to LRS-Gowdy 



The Gowdy T 3 model with linear polarization is symmetric under the interchange of the 
directions coordinatized by a and 5. Owing to this, it has a subset of classical solutions 
with local rotational symmetry (LRS), in which the scale factors of these two directions can 
be identified during the entire evolution. We can then restrict the Gowdy model, both in 
vacuo and with matter, to the LRS-Gowdy model in which every solution is of this kind. The 
restriction can be performed classically, prior to quantization, or starting with the quantized 
model. We will focus our attention on the latter approach, passing from quantum Gowdy to 
quantum LRS-Gowdy, and leave for the interested reader the proof that the quantum model 
obtained in this way is indeed recovered by a direct quantization of the classical LRS-Gowdy 
spacetimes along the lines explained in this work. 

In analogy with the discussion of Ref . [l6| , in which the quantum FRW model is obtained 
from quantum Bianchi I, we define the following map from (generalized) states associated 
with the Gowdy model to those of the LRS-Gowdy cosmologies, denoted by \if)(v,A)): 



|%A,T)) 



(3.15) 



T 



The sum is carried out over all values of T in the considered superselection sector. Applying 
this map in the Hamiltonian constraint ( I3.10p . we obtain 



4 

7fG 



d 



\^(v,A))+ £ £ a£(*M + A)) 

kS{0,4} se{+,-} 

SK + V 



+ 



2A 



fte{o,4} se{+,-} 
Hl + H$] |V>(v,A)> = 0, 



x s K {v)\ip (sk + v, A + w v (sk))} 



(3.16) 



where \ip s K (sK + v, A)) are the combinations 

\ijj s K (sk + v, A)) = \if> (sk + v, A + w v (sk) - w v (s2))) 

+ \ip {sk + v, A + w v {sk))) + |V> {sk + v, A + w v {s2))). 



(3.17) 



As we have already remarked, the result agrees with the constraint obtained by a suitable 
hybrid quantization of the classical LRS-Gowdy model. It is worth noting that the intro- 
duced map works because the coefficients appearing in Eq. (I3.10p do not depend on the 
variable T, over which one sums to perform the projection. Indeed, this kind of map only 
makes sense if the classical model admits the imposition of an additional symmetry which 
allows its reduction into a dynamically stable submodel. A similar projection summing over 
A is not viable because the coefficients of the inhomogeneous contributions in the constraint 
depend explicitly on this variable, reflecting the fact that the associated kind of isotropic 
solutions exist just when the inhomogeneities are unplugged. 



IV. DISCUSSION 



We have completely quantized the Gowdy T 3 model with linearly polarized gravitational 
waves provided with a minimally coupled massless scalar field as matter content. The 
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description adopted for the matter inhomogeneities is such that they can be treated in 
exactly the same way as the gravitational ones, the former just duplicating the contributions 
of the latter in the constraints. In this situation, we have been able to apply the hybrid 



quantization methods developed in Refs. [191-42 21] almost straightforwardly to this system 



with local physical degrees of freedom both in the matter content and in the gravitational 
field. To our knowledge, it is the first time that a model with these properties has been 
quantized to completion in the framework of LQC 

Since the structure of the constraints when the matter field is present is the same as 



in vacuo, all the results obtained in Ref. [21] for the vacuum Gowdy model apply as well 



to our model. Thus, in particular, we recover on physical states the standard quantum 
field theory description of both the matter and the gravitational inhomogeneities, living 
on a cosmological background quantized using LQC methods and consisting of a Bianchi I 
universe with a homogeneous massless scalar field. In addition, it is guaranteed that the 
states which are the analog of the classical singularity decouple naturally in the quantum 
model, so that, to this extent, the initial singularity is resolved at the kinematical level. 

Conceptually, the hybrid quantization of the present family of inhomogeneous cosmologies 
has introduced no technical complication with respect to the vacuum case. Nonetheless, the 
situation is radically different when one considers the interest of the quantum model from 
a physical point of view. In fact, thanks to the inclusion of the massless scalar field, the 
homogeneous sector of the Gowdy T 3 model, namely the Bianchi I model, admits now 
isotropic flat solutions of FRW type, while in vacuo only the trivial Minkowskian solution 
is allowed. 

On the other hand, the analysis of the classical solutions of the linearly polarized Gowdy 
T 3 model in vacuo [35] and the study of the effective dynamics obtained from the hybrid 
quantization of this model [36] show that small inhomogeneities do not increase arbitrarily 
in the evolution. Then, if we consider initial data which are sufficiently close to homogeneity, 
the corresponding solution would remain approximately homogeneous during the evolution. 
Besides, in the nonvacuum model, matter and gravitational inhomogeneities evolve in iden- 
tical way. This strongly indicates that initial data in a sufficiently small neighborhood of 
those with isotropy and homogeneity have to lead to approximately isotropic and homoge- 
neous solutions. Therefore, it is natural to compare the dynamics of our Gowdy model with 
that of the flat FRW model (with three-torus topology) in order to see how the inclusion of 
anisotropies and inhomogeneities affects the evolution of a flat FRW background. Moreover, 
we are now in a perfect situation to carry out this comparison at the quantum level, since 
here we have accomplished the full quantization of the Gowdy T 3 model in the presence of 
the massless matter, and the loop quantization of the FRW model coupled to the homoge- 
neous massless field is well known p, 0] . Even though the inhomogeneities in our model are 
not all those allowed in a universe like the one which we observe (but just a subfamily with 
the symmetries of the Gowdy T 3 cosmologies), their analysis should shed light on the kind 
of quantum effects affecting the evolution and on the consequences of the quantum geometry 
on the primordial fluctuations. 

For these purposes, it is preferable to focus on the LRS-Gowdy model derived in Sec. 
IIII CI Indeed, the consideration of the two degrees of freedom of anisotropy that the homo- 
geneous sector of the general Gowdy model possesses would only complicate the equations 
unnecessarily. The presence of either two degrees or just one degree of anisotropy does not 
seem to have any conceptual relevance for the proposed analysis. 

In order to face this analysis, the idea is to add and subtract in Eq. (13.161) the term 
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corresponding to the FRW model, which coincides with the first line of Eq. ( 13 . 16[) but keeping 
the variable A unchanged. We can then rewrite Eq. (13.161) as the constraint equation of the 
FRW model coupled to a homogeneous massless scalar field plus a number of contributions 
coming from all other terms. These contributions contain the inhomogeneities and the 
difference between the gravitational parts of the constraints for LRS-Bianchi I and FRW, 
a difference which is due to the anisotropics. In this way, the resulting expression modifies 
the densitized Hamiltonian constraint of the FRW model by the effects of the anisotropies 
and inhomogeneities, so that it is not longer equal to zero. As we have commented, we 
are interested in comparing the FRW model with the inhomogeneous LRS-Gowdy model 
when these inhomogeneities and anisotropies are small. In this regime, it makes sense to 
apply a type of Born-Oppenheimer approximation, similar to others commonly employed 



in cosmology (see e.g. [371 138|] ) , and assume that the variations of the isotropic degrees of 



freedom and those of the rest of degrees have considerably different typical scales, therefore 
giving them a different status. Then, in this approximation, it is easier to derive effectively 
the influence that anisotropies and inhomogeneities produce on the isotropic background. 
We leave for future research this detailed analysis. 



ACKNOWLEDGEMENTS 

The authors are grateful to L. Garay and J. Olmedo for discussions. This work was 
supported by the Spanish MICINN Project No. FIS2008-06078-C03-03 and the Consolider- 
Ingenio Program CPAN No. CSD2007-00042. M. M-B and D. M-dB have been supported by 
CSIC and the European Social Fund under the grants I3P-BPD2006 and JAEPre_09_01796, 
respectively. 



[1] M. Bojowald, Living Rev. Relativity 11, 4 (2008). 

[2] A. Ashtekar, Nuovo Cimento B122, 135 (2007). 

[3] G.A. Mena Marugan, AIP Conf. Proc. 1130, 89 (2009). 

[4] A. Ashtekar, M. Bojowald, and J. Lewandowski, Adv. Theo. Math. Phys. 7, 233 (2003). 

[5] A. Ashtekar and J. Lewandowski, Classical Quantum Gravity 21, R53 (2004). 

[6] T. Thiemann, Introduction to modern canonical quantum general relativity (Cambridge Uni- 
versity Press, Cambridge, England, 2007). 

[7] A. Ashtekar, T. Pawlowski, and P. Singh, Phys. Rev. D 74, 084003 (2006). 

[8] A. Ashtekar, A. Corichi, and P. Singh, Phys. Rev. D 77, 024046 (2008). 

[9] M. Martin-Benito, G.A. Mena Marugan, and J. Olmedo, Phys. Rev. D 80, 104015 (2009). 

[10] A. Ashtekar, T. Pawlowski, P. Singh, and K. Vandersloot, Phys. Rev. D 75, 024035 (2007). 

[11] L. Szulc, W. Kaminski, and J. Lewandowski, Classical Quantum Gravity 24, 2621 (2007). 

[12] K. Vandersloot, Phys. Rev. D 75, 023523 (2007). 

[13] E. Bentivegna and T. Pawlowski, Phys. Rev. D 77, 124025 (2008). 

[14] D. Chiou, Phys. Rev. D 75, 024029 (2007). 

[15] M. Martm-Benito, G.A. Mena Marugan, and T. Pawlowski, Phys. Rev. D 78, 064008 (2008). 

[16] A. Ashtekar and E. Wilson-Ewing, Phys. Rev. D 79, 083535 (2009). 

[17] A. Ashtekar and E. Wilson-Ewing, Phys. Rev. D 80, 123532 (2009). 

[18] E. Wilson-Ewing, Phys. Rev. D 82, 043508 (2010). 



12 



[19] M. Martin-Benito, L.J. Garay, and G.A. Mena Marugan, Phys. Rev. D 78, 083516 (2008). 
[20] G.A. Mena Marugan and M. Martin-Benito, Int. J. Mod. Phys. A 24, 2820 (2009). 
[21] L.J. Garay, M. Martin-Benito, and G.A. Mena Marugan, Phys. Rev. D 82, 044048 (2010). 
[22] M. Martin-Benito, G.A. Mena Marugan, and E. Wilson-Ewing, Phys. Rev. D 82, 084012 
(2010). 

[23] R.H. Gowdy, Ann. Phys. 83, 203 (1974). 

[24] M. Bojowald, Classical Quantum Gravity 20, 2595 (2003). 

[25] A. Corichi, J. Cortez, and G.A. Mena Marugan, Phys. Rev. D 73, 041502 (2006). 
[26] A. Corichi, J. Cortez, and G.A. Mena Marugan, Phys. Rev. D 73, 084020 (2006). 
[27] A. Corichi, J. Cortez, G.A. Mena Marugan, and J.M. Velhinho, Phys. Rev. D 76, 124031 
(2007). 

[28] A. Corichi, J. Cortez, G.A. Mena Marugan, and J.M. Velhinho, Classical Quantum Gravity 
23, 6301 (2006). 

[29] J. Cortez, G.A. Mena Marugan, and J.M. Velhinho, Phys. Rev. D 75, 084027 (2007). 
[30] J.F. Barbero G., I. Garay, and E.J.S. Villasehor, Phys. Rev. D 74, 044004 (2006). 
[31] J.F. Barbero G., D. Gomez Vergel, and E.J.S. Villasehor, Classical Quantum Gravity 24, 5945 
(2007). 

[32] M. Pierri, Int. J. Mod. Phys. D 11, 135 (2002). 

[33] A. Corichi, J. Cortez, and H. Quevedo, Int. J. Mod. Phys. D 11, 1451 (2002). 
[34] J. Cortez and G.A. Mena Marugan, Phys. Rev. D 72, 064020 (2005). 
[35] V. Moncrief, Phys. Rev. D 23, 312 (1981). 

[36] D. Brizuela, G.A. Mena Marugan, and T. Pawlowski, Classical Quantum Gravity 27, 52001 
(2010). 

[37] J.J. Halliwell and S.W. Hawking, Phys. Rev. D 31, 1777 (1985). 

[38] C. Rovelli and F. Vidotto, Classical Quantum Gravity 25, 225024 (2008). 



